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1 Introduction 



In this note we discuss transverse momentum distributions in heavy flavour decays. As may be easily recognized, 
the structure of the expressions reflects the analogous case of the energy-energy correlations or shape variables 
in e + e~ annihilation, i.e. it involves a coefficient function, a universal function / containing the logarithmic 
contributions and a remainder function. This structure shows, on the other hand, the general form common to 
all the reactions such as Drell-Yan (DY), Deep Inelastic Scattering (DIS), Jet Shape variables, etc., once the 
large logarithmic contributions are separately resummed in the function /. 

The purpose of computing the set of perturbative contributions for the heavy quark decay case is introduced 
in this paper together with some preliminary results. We compute the next-to-leading order (NLO) coefficient 
B\ describing single-logarithmic effects, which, together with the universal A\ and A2 coefficients, determines 
the function / to NLO accuracy. 

We outline here the peculiar phenomenological content of the transverse momentum distribution as making, 
for several reasons, a stand-alone case in heavy flavour physics. 

The main phenomenological application of our work involves the spectrum of tranverse momenta of the 
produced hadrons, with respect to the photon direction, in the reaction b — > sj. This process appears to be 
very clean because, in the rest frame of the B-meson, the photon momentum unambigously fixes the reference 
direction for transverse hadron momenta. 

The final goal of our study is a comparison of the complete NLO perturbative distribution with equally 
accurate experimental data. Such a comparison should allow the extraction of the non-perturbative component 
in the process i.e. the Fermi motion effects of the &-quark inside the B-meson and the final-state hadronization 
of the strange quark. The information about non-perturbative dynamics obtained with p± distributions is 
complementary to the one obtained with the well-studied threshold distributions. Roughly speaking, transverse 
momentum distributions are sensitive to the motion of the heavy quark in the plane orthogonal to the decay 
axis, while threshold distributions are sensitive to the motion along the decay axis. This topic deserves some 
attention since, as well known, the transverse momentum distributions of heavy flavours in hadronic processes 
are presently poorly understood. 

The paper is organized as follows. In sec. 2 we present our result on transverse momentum distributions. 
In sec. 3 we summarize the main results on threshold distributions. In sec. 4 we discuss the results and we 
compare the two different distributions. Finally, sec. 5 contains the conclusions together with an outlook at 
future developments. 



2 Transverse momentum distribution 



As anticipated in the introduction, let us consider the distribution of relative transverse momenta between a 
strange hadron h s and the photon in radiative B decays: 

B^h s +X + j, (1) 

In practice, we expect h s to be a meson such as a K, a K* , etc.. Without any generality loss, we can take 
the B-meson at rest and the photon flying along the plus direction (+z axis); the relative tranverse momentum 
then coincides with the projection of the h s momentum in the x — y plane. In general, we may identify three 
different mechanisms in transverse momentum dynamics: 

1. soft-gluon emission. The elementary process in (|l|) is: 

b s + X + 7. (2) 
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Non-vanishing tranverse momenta 

P v + (3) 

are generated by soft-gluon emissionR 

2. Fermi-motion. The beauty quark and the light valence quark in the £?-meson exchange soft momenta with 
each other. This implies that, before decay (Q) occurs, the b quark has a non-zero transverse momentum 

P/ ~ A (4) 

where A is the QCD scale. 

3. final-state hadronization. The strange quark emitted in the hard process fragments into the final hadron 

h s , 

s — ► h s + (anything) . (5) 

Since the latter is a soft process, this implies a change in the transverse momentum of the strange system 
by a quantity 

Ph ~ A. (6) 

The total tranverse momentum of h s is then 

P = P P +P/ + Pft- (7) 

Let us observe that Fermi motion and hadronization effects are of the same order and cannot therefore be 
separated from each other. As we are going to show later, the situation is instead different in threshold 
distributions. According to the above estimates, if we take 

|p| » A, (8) 
we have that the non-perturbative effects described in 2. and 3. are negligible, so that 

P ^ P P , (9) 

and a (resummed) perturbative computation controls the distribution. The transverse momentum distribution 
in impact parameter space b jl], has the characteristic form || j|: 

^^r=C(a s ) f(b;a s )+R(b;a s ). (10) 
1 b do 



Tb is the Born width and is given by: 



7T 327T J 

where we consider only the leading operator O7 in the effective b — > sj hamiltonian . 

The coefficient function C(as) is process-dependent and short-distance dominated and has an expansion in 
powers of as of the form: 

C(a s ) = l + ^c+(^) 2 c' + ... (12) 



1 Transverse momenta are generally denoted with boldface symbols 
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where Cp — (N% — l) / (2N C ) = 4/3, as = as{Q 2 ) and Q = m B is the hard scale. The remainder function R 
is process-dependent and has an expansion starting at order as of the form 



R(b;as) = ^r(b) + (^) 2 r>(t) + .... 



\ 7T 

R vanishes in the limit of a large impact parameter: 

R(b;as)^0 for b — » oo, 
implying that it is short-distance dominated. 



(13) 



(14) 



The function / (6; as) contains the large logarithmic contributions in a resummed form and can be written 
as the exponential of a series of functions JO: 



/ (b;a s ) = exp [L#i(/3 asL) + g 2 {fioasL) + a s g^PoasL) + 



where 



-';,2 



L = lot 



68 



(15) 



(16) 



and bo = 2 exp [— 7_e] — 1.12 with ^p — 0.577 the Euler constant. Our computation of the functions g\ and g 2 
gives: 



A\ 1 

.91 M = ^ — [log(l-w)+a;] 
2p w 



1 - w 



+ log(l-w) 



2/3 3 



lQg(l - LU) 
1-LU 



W 1,9/ 

+ -log 2 

1 — u> 2 



(17) 



B 



-=ilog(l- W ),(18) 
Po 



where 



7T 



and 



A 2 = ^ 



Ca 1 36 ~ 12 



5 t ' 



(19) 



The above value of A 2 is in the MS scheme for the coupling constant j|. As usual, Ca — N c = 3, Tr = 1/2 
and rip = 3 is the number of active quark flavours. The first two coefficients of the /3-function are: 



(3o = 



llC A -2n F 33-2n F 



12tt 



12tt 



Pi = 



YIC\ - 5C A n F - 3C F n F _ 153 - 19 n F 



24tt 2 



24tt 2 



(20) 



A complete NLO analysis requires also the knowledge of the coefficient function c and of the remainder function 
r (b) of order as, which at present are unknown. Their computation is in progress and will be presented in a 
forthcoming publication. 



Denoting by Si the soft contribution at one loop and by C\ the collinear one in the notation of ref. 



ft-^, 



r 3C F 



the next-to-leading coefficient B\ is given by 



Bi = 



Si 



c x = - 



5CV 

4 7T 



(21) 



(22) 



In usual hard processes, such as DIS or DY, Si = so that B\ — C\. In our case, single logarithmic effects are 
more pronounced because this coefficient is almost a factor of 2 larger. 
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The expansion to order a 2 s of the exponent reads: 



1 



log/ (6; a s ) = -^A ias L 2 



(23) 



B x a s L - i^AjalL 3 - jA 2 a 2 s L 2 - X -B^a 2 s L 2 . 

Let us note that a single constant B\ controls the single-logarithmic effects in any order. The physical reason 
is that a soft gluon and a collinear one with the same transverse momenta are emitted with the same effective 
coupling as (kj) This is to be contrasted with the threshold case (see next section). The function f(b, as) 
is plotted in fig. 1. 



Fig. 1: Plot of the junction f(b) in the variable y — Q b jb\ (as — 0.22). Solid line: NLO; dotted line: LO. 



The resummed distribution ( |l0| ) becomes singular when 

ui -f r 

Since 

lo = p a s L 



logQ 2 6 2 



(24) 



(25) 



logQ 2 /A 2 ' 

the singularity occurs when the transverse strange momentum becomes as small as the hadronic scale, 

p± w i « A. (26) 

The singularity ( p4| ) is produced by the infrared pole in the running coupling and signals an intrinsic limitation 
of resummed perturbation theory in agreement with previous qualitative analysis. Let us note that the function 
g 2 has basically a pole singularity in the limit (E4j), while <?i has only a softer, logarithmic, singularity. 



3 Threshold distribution 



In this section we recall the main results on threshold distributions. The distribution in Merlin space — the 
TV-moment of the rate — has a similar representation to ( fl0| ) 

-Lr N =C(a s ) In (as) + Rn (a s ) , (27) 

1 B 
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where now the large logarithm contains the JV-variable: 

N 

L = log — (threshold case) , (28) 

N 

with No = exp [—7b] . The functions <?i in the exponent are different with respect to the ones in the p± case 
and the leading and next-to-leading ones read jfO], [llj : 

9i (A) = -^[(l-2A)log(l-2A)-2(l-A)log(l-A)], 



92 (A) = P ° A \-f lAl [log(l - 2A) - 21og(l - A)] - ^ [log 2 (l - 2A) - 21og 2 (l - A)] 



+Al og( i_2A) + ^log(l-A). (29) 

The expansion to order a s of the exponent reads: 

log / w = - X -A x a s L 2 - a s (^ + C\) L — ^Aifl Q a s L 3 - ^A 2 a 2 s L 2 - (sj. + ^C^j l3 a 2 L 2 . (30) 
Let us comment on the above results. The single-logarithmic effects at one loop are controlled by the constant 

S l + C 1 = ~% (31) 

4 7T 

i.e. by the sum of the soft and the collinear coefficients, which is different from the p^-case (cf. eqs. (|22|) and 
(|2^)). At two- loop they are instead controlled by a different constant, 

z on 

The soft and the collinear terms begin to differentiate at this order and the soft one has a two times larger 
coefficient. Contrary to the p± case, two different constants are needed to describe the single logarithmic effects. 
The dynamical difference between soft and collinear terms is that, for a fixed jet mass, the transverse momentum 
of a soft gluon is substantially smaller that of a collinear gluon P, [121 pi pT[ . 



The functions g\ and g 2 in (|29J) — and therefore also the resummed distribution — have two different 
singularities 0, || § : 

i) the first one occurs when 

1 _ log Q 2 /m 2 

2 ~ logQVA^' [d6) 

or, equivalently, when 

m 2 w AQ, (34) 

where m is the mass of the final hadronic jet s + X. In the last member of (|3^), we have used the 
approximation TV « Q 2 /m 2 [ ^6| . The singularity (j33|) signals the occurence of non-perturbative effects in 
region ( p4| ) — to be identified with the well-known Fermi motion jlTj]; it is related to soft-gluon effects 
— i.e. to the terms proportional to At, Si and A 2 — and not to collinear ones — the term proportional 
to C\. Fermi- motion effects are therefore controlled by soft and not by collinear dynamics. This fact 
allows a factorisation of Fermi-motion effects by means of a function taking into account soft dynamics 
only, the well-known shape functior^j |Q . In this region, initial bound state effects become relevant while 
final-state binding effects can be neglected Jl9[ |[ pj . 



The shape function is also called structure function of the heavy flavours. 
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ii) the second singularity occurs at 



A = l, 



(35) 



or 



m' 



2 



A 2 



(36) 



and is related to final-state hadronization effects. Both soft and collinear terms are singular in this region 
and there are non-perturbative effects related to initial as well as final bound-state dynamics. 



4 Discussion and comparison of the distributions 

Let us now discuss the physics of p± distributions and compare with the threshold case. In the latter case, 
there is a first singularity, of soft nature, closer to the origin in A = 1/2, which can be removed by introducing 
the shape function. In the p± case, an analogous singularity is absent, as also supported by physical intuition. 
In both distributions, there is a singularity in one (A = 1 or u> = 1), which cannot be removed by introducing 
a non-perturbative soft function. This implis that, in the p± case, non-perturbative effects cannot be treated 
with a "shape-function" approach. One could try a different approach based on the effective theory introduced 
in |l3| or the collinear effective theory developed in [ po| . 

Let us note that, in general, the singularities of the functions gi (u>) are more severe than those of the functions 
Qi (A) for A — > 1/2 or 1. For example, g\ (ui) has a logarithmic singularity, while g\ (A) has an additional prefactor 
1 — 2A or 1 — A which softens the singularity. Owing to the different singularity structure, the pj_-distribution 
is complementary to the threshold one and gives independent information about non-perturbative physics. 

Finally, let us make a general field-theory remark. Threshold distributions and p^-distributions have a 
different theoretical status. The former are completely inclusive quantities and can be computed as the imaginary 
part of a forward scattering amplitude: a non-perturbative computation with lattice QCD is in principle feasible 
pl| . This is not the case for p^-distributions, which instead are true jet quantities: a separate computation of 
real and virtual diagrams is unavoidable and a lattice QCD computation is in principle unfeasable. 



5 Conclusions 

In this note we have discussed transverse momentum distributions in radiative -B-decays and we have computed 
the next-to-leading coefficient B\ . The structure of the logarithmic corrections is analogous to the one encoun- 
tered in energy-energy correlations or shape variable distributions in e + e~ annihilations. Owing to the different 
singularity structure with respect to threshold distributions, it does not seem possible to define the analogue of 
a shape function. An operator definition of the non-perturbative effects in this case should presumably involve 
a different effective theory explicitly containing the transverse degrees of freedom. 

The comparison of our NLO distribution with accurate experimental data may give new and independent 
information about the effective size of the non-perturbative corrections of order A/ms; a recurrent problem 
in B-physics is indeed the separation of perturbative effects from non-perturbative ones and the estimate of 
the latter. The computation of the remaining next-to-leading terms — the one loop coefficient function and 
remainder function — is in progress and will be presented elsewhere. 
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